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Abstract 

0^ ' The U-CAS is a spinning magnetized top that is levitated in a static magnetic field. The field is produced by a 
^^ , permanent magnet base, positioned below the hovering top. In this paper we derive upper and lower bounds for ft,m~the 
maximum hovering height of this top. We show that the bounds are of the form alo where a is a dimcnsionless number 
ranging from about 1 to 12, depending on the constraints on the shape of the base and on stability considerations, and 
Iq is a characteristic length, given by 
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Here, /io is the permeability of the vacuum, fi is the magnetic moment of the top, Mq is the maximum magnetization 
of the base, m is the mass of the top, and g is the free-fall acceleration. For modern permanent magnets we find that 
Iq ~ 1 [meter], thus limiting h to about few meters. 

Index Terms- U-CAS, Levitron, magnetic trap, magnetic levitation, hovering magnetic top. 

I. Introduction 
A. What IS the U-CAS? 

The U-CAS is an ingenious device that hovers in mid-air while spinning. It is marketed as a kit 
in Japan under the trade name U-CAS 1^, and in the U.S.A. and Europe under the trade name 
Levitron"""^ |0], |^, [Q. The whole kit consists of three main parts: A magnetized top which weighs 
about 18 gr, a thin (lifting) plastic plate and a magnetized square base plate (base). To operate the 
top one should set it spinning on the plastic plate that covers the base. The plastic plate is then raised 
slowly with the top until a point is reached in which the top leaves the plate and spins in mid-air above 
the base for about 2 min. The hovering height of the top is approximately 3 cm above the surface of 
the base whose dimensions are about 10 cm x 10 cm x 2 cm. The kit comes with extra brass and 
plastic fine tuning weights, as the apparatus is very sensitive to the weight of the top. It also comes 
Q with two wedges to balance the base horizontally. 

J>-^ B. The adiabatic approximation. 

^ '. 
p i. The physical principle underlying the operation of the U-CAS relies on the so-called 'adiabatic 

^ I approximation' f^, |^, 0: As the top is launched, its magnetic moment points antiparallel to the 

"^ ' magnetization of the base in order to supply the repulsive magnetic force which will act against the 

'^ gravitational pull. As the top hovers, it experiences lateral oscillations which are slow ( f^iaterai — 1 

5^ ,Hz) compared to its precession (f^precession ~ 5 Hz). The latter itself, is small compared to the top's 

spin (figpin ~ 25 Hz). Since f^spin ^ ^precession the top is considered 'fast' and acts like a classical 

spin. Furthermore, as f^precession ^ ^lateral this spin may be considered as experiencing a slowly 

rotating magnetic field. Under these circumstances the spin precesses around the local direction of 

the magnetic field B (adiabatic approximation) and, on the average, its magnetic moment /x points 

antiparallel to the local magnetic field lines. In view of this discussion, the magnetic interaction energy 
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which is normally given by —fi ■ B is now given approximately by /i|B|. Thus, the overall effective 
energy 'seen' by the top is 

Ees{r) c:^mgz + fi\B{r)\, (1) 

where m is the mass of the top, g is the free-fall acceleration and z is the height of the top above 
the base. By virtue of the adiabatic approximation, two of the three rotational degrees of freedom 
are coupled to the transverse translational degrees of freedom, and as a result the rotation of the axis 
of the top is already incorporated in Eq.(|l|). Thus, under the adiabatic approximation, the top may 
be considered as a point-like particle whose only degrees of freedom are translational. The important 
point of this discussion is the following: The energy expression written above possesses a minimum 
for certain values of [i/mg. Thus, when the mass is properly tuned, the apparatus acts as a trap, 
and stable hovering becomes possible. A detailed description of this device, extending beyond the 
adiabatic approximation, may be found elsewhere |0, 0, p. For the purpose of this paper, the 
adiabatic approximation will suffice. 

C. The limitations on the hovering height. 
In this paper we focus on the question: 

How high can the hovering height of the U-CAS be? 

First, we define what do we mean by 'height': We assume that the top hovers above some horizontal 
plane. The height of the top is measured with respect to that plane. The 'rules of the game' are, that 
we can put permanent magnet helow the plane, but never above this plane. Whatever is below this 
plane, will be called 'the base'. The method we use to answer the question above will be by getting 
upper and lower hounds for this height, denoted by km- We begin by pointing out the factors that 
limit the hovering height of the U-CAS. 

First, we set aside the question of stability and assume that the top is guided by a vertical axis, 
allowing it to move only axially (up and down). In this case, the easiest way to increase /i-the hovering 
height, is by using a more powerful magnet for the base. This technique however, cannot be applied 
indefinitely since there is no way to increase the magnetization of a given substance without limit. 
We are therefore forced to limit the strength of M(r)-the magnetization density of the base, to some 
maximum value, say Mq, and design the shape and magnetization density of the base so as to maximize 
h. One more factor that limits the levitation height is the amount of magnetic substance, or volume, in 
our disposal. Clearly, the more material we use, the larger is the levitation height that can be achieved. 
But still, even if the volume is infinite the levitation height is bounded, since |M(r)| is bounded. 

When stability is taken into account, things get more complicated. Stability sets another limitation 
on the design of the base, in addition to the limitations that were discussed above: A brief look at 
Eq.(P tells us that for the top to hover stably over the plate, the effective energy Ecs should possess 
a minimum. This means, in particular, that 



d^E,s{r) 52|B(r)| d''\B{r)\ 9^ |B (r) 



dv'^ dx^ dy^ dz 



> 0, (2) 



where the derivatives are evaluated at the equilibrium position of the top. Since these are homogenous 
inequalities, it is clear that the region in space where a minimum may occur, does not depend on the 
strength of the magnetic field. As a consequence, the hovering height of the U-CAS is not determined 
by the strength of the magnetic field but on its geometry, or alternatively, by the shape of the base. 



As an example, it has been shown |T^ that for a vertically magnetized base in the shape of a disk of 
radius R, the range of heights h for which stable hovering is possible, is very narrow and is around 
h ~ R/2. This result agrees roughly with the parameters of the U-CAS, for which /? ~ 10 cm and 



h ^ 4: cm, when h is measured from the center-of-mass of the base. The strength of the field then 
comes into play by tuning it in order to achieve equilibrium for a given mass of the top. Equilibrium 
prevails when the total force on the top vanishes, i.e. when 

F= -V^eff = -"^c/z-/iV|B(r)| = 0. (3) 

Both in the guided case and in the stable case, we see that in order to increase the hovering height 
of the U-CAS, we should design a better base. We do expect however, that the hovering height cannot 
be increased indefinitely. In connection with the redesigning of the base, we have recently shown [TT|, 
both theoretically and experimentally, that the use of a vertically magnetized ring of radius i? as a 
base, increases the hovering height by more than three times to about 1.7i?!. This increase in height 
did not came without cost, as the tolerance on the mass of the top became more stringent, being 



Am/m ~ 0.6% for the ring vs. Am/m ~ 1.2% in the case of the disk flOl . 

D. The structure of this paper. 

All the calculations that we do are outlined in Sec.(|I|). We will now describe what we calculate and 
what is the motivation behind it. 

We start in Sec. ( [11- A] ), by deriving a close form expression for the first and second derivatives of the 
magnetic field, in terms of the magnetization density of the base. The first derivative is the magnetic 
force on the top, which is used throughout the paper, while the second derivative is exploited when 
we discuss stability. 

In Sec.( [II-BD we consider the problem of maximizing the hovering height in the case where the volume 
of the base may be infinite, and may be arbitrarily magnetized under the constraint that |M(r)| = Mq. 
We still do not require that the top be stable against lateral translations, and assume that it is guided 
along a vertical axis. We show that in this case 

'''max — J-'^'O) 

where /q is a characteristic height for the problems that are discussed in this paper. It is defined as 

_ /io /iMo 
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where /ig is the magnetic permeability of the vacuum, /i is the magnetic moment of the top, m is the 
mass of the top, and g is the free-fall acceleration. 

In Sec.( [lI-CD we consider the same problem as before but for a base which is uniformly magnetized 
along the vertical direction, namely M(r) = Mqz, where z is a unit vector in the vertical direction. In 
this case we find that 

'^max — O.OtQ- 

In order to arrive to better bounds on km, our next step is to limit the amount of material available 
for the base. Thus, in Sec.( P-D| ) we maximize the hovering height in the case where the base may be 
arbitrarily magnetized under the constraint that |M(r)| = Mq, and that its volume Vq is given. Here, 
the result is given in the form of a plot, showing the dependence /imax on Vq- In the limit Vq —>■ oo. 



we recover the result of Sec. ([II-B|) . Our derivation is also constructive in the sense that it shows how 



to construct the optimal base for a given Vq. In this paper however, we do not discuss this matter 
thoroughly due to space shortage. 

For completeness, we study in Sec. (|II-ED the same problem as in Sec.( |II-D]) , but for a uniformly 
magnetized base. Here, again we find the dependence of /imax on Vq and recover the result of Sec.(|I] 

g)- 



Sec.( II-F ) is the first part where stabihty is added into play: We begin by explaining how to test for 



stability of the top under the adiabatic approximation, and study the possible levitation heights of a 
uniformly magnetized base whose shape is cylindrical of volume Vq. As the stability condition is given 
in the form of inequality relations (see Eq.(|^)), it is not suffice to determine uniquely the levitation 
height. We therefore choose to study two particular stable points: an isotropic stable point and the 
highest possible stable point. The meaning of these points would become clear later. In this case we 
also give a plot, showing the dependence of /imax on Vq. This result should also be considered as a 
lower bound for hm, since the cylindrical base is a special case of all the base configurations that are 
possible. 

In section( [ri-G|) we take one step further, and generalize the result of section (p.I-I']) to the case of a 



base which may be arbitrarily magnetized, and look for the lowest possible stable point. In this case 
however, we assume that Vq = cxd, to ease the solution of the problem. We also show that in this case, 
where /imax — 8.4/o and Vq = oo, the lowest possible stable point and the highest one coincide, so this 
result represent not a bound for hm, but is actually hm itself, for the case Vq = cxo. 

In Sec. ( |IIip we discuss interesting aspects of our results. In particular we estimate the value of /q 



for modern permanent magnet materials and calculate the values of the bounds for hm- We comment 
on the implications of our results and discuss other related questions. 

II. Mathematical Formulation 
A. The upward magnetic force and its derivative. 

A simplified model of the U-CAS is shown in Fig.(|l|). It consists of a point-like particle of mass 
m and magnetic moment fi (pointing downward), hovering at a height h above the z = plane. For 
the moment, we consider the top as if it was guided along the z-axis so that only its vertical motion 
is allowed. Other degrees of freedom are considered 'frozen'. The region z < may be partially or 
wholly filled with a magnetized substance ("base"), whose magnetization density is denoted by M(r), 
producing a magnetic field throughout space which we denote by B(r). In what follows, we assume 
that the base (and hence the magnetic field) has a cylindrical symmetry around the z-axis. Thus, along 
this axis the magnetic field possesses only a z-component. We further assume that this component is 
positive, namely that the magnetic field along the 2;-axis is pointing upward. 

Using Eqs.(|l]),(|D we find that, when the top is in equilibrium, the total vertical component of the 
force on the top vanishes, 

5|BL 
Fz = -mgz - yu^— z = 0, 
oz 

where the term —mgz is the gravitational force pulling the top towards the base, and the term 
~^zd\^\ /dz is the magnetic force which pushes the top upward. Note that d\^\/dz should be 
negative. 

Since the top is allowed to move only along the 2;-axis, we may write |B| = Bz (in what follows 
we deserve the notation Bz to denote the 2;-component of the magnetic field along the 2;-axis, namely 
Bz = z ■'B {x = 0,y = 0,z)). Now, the magnetic force on the top simplifies to 

J, dBz. 

oz 



and hence, when the top is in equilibrium aX z = h 

OB, 



-/i 



dz 



mg. (4) 



z=h 

Reciprocity allows us to express dBz/dz in terms of M(r) as 



dz 



Mfr') 



dB„ 



^) ^3 / 



dz 



d'r'. 



(5) 



Here, Brf(r'; z) is the field at the point r' produced by a unit magnitude dipole pointing upward, located 
at a height z along the 2;-axis. Taking r and 6 as depicted in Fig.(|l|), we can write B^ as 



B. 



/io (3cos6'f — z) 
47r r^ 



(6) 



where f is a polar unit vector, also defined in Fig.(|ip, and /xq is the magnetic permeability of the 
vacuum. It is worth to note that d^d/dz is nothing but the field produced at r' by a quadruple 
located at a height z along the z-axSs. This quadruple is made out of a pair of identical dipoles: One 
is located at z and pointing downward and the other is at an infinitesimally higher position z + dz and 
pointing upward, with their magnetic moment being 1/dz. 

When Eq.dl) is substituted into Eq.(|) we find that, dX z = h, the field derivative is given by 



dB.. 



dz 



^ I sin 6de 

47r 



1 



27rr^dr—M{r,e) •n(0). 



(7) 



7r/2 



-h/cos{e) 



where 



n(^) = 6 sin ^ cos 99 + 3(3 cos^ 9 - l)f , 



and where ^ is a polar unit vector orthogonal to f, as is shown in Fig.(|l]). 

We would also need the second derivative of the field, when we discuss stability. It is given by 



d'B., 



dz^ 



Ait 



sin 9d9 



2'Kr'^dr—M{r,9) ■s{9), 



(9) 



7r/2 



-h/cos{e) 



where 



s{9) = 9 sin ^ (5 cos^ 9-1)9 + 12 cos ^(5 cos^ 9 - 3)f . 



(10) 



Note that in both Eqs.(0) and (|), M(r, 9) is defined with respect to a polar coordinate system whose 
origin is at z = h and not at 2; = 0. Also note, that in going from Eq.(|^ to Eq.(^, one should take 
into account the dependence of the unit vectors 9 and f on z. 

B. Arbitrarily magnetized, infinite base. 

We consider an infinitesimal volume element at some point (r, 9) within the base. It is essentially 
a magnetic dipole whose magnitude is Mq, and whose direction we will now find by the requirement 
that the magnetic force on the top is maximized. We have already shown that the force on the top, 
contributed by one elemental dipole, is proportional to the magnetic field B^^ that would be produced 
at (r, 9) by a quadruple located at z = h. To make this force maximal we assign to each point (r, 9) 
in 2 < 0, a magnetization density M(r) whose magnitude is Mq and whose direction is antiparallel to 
the field line of B^ at that point. The magnetization density so defined will maximize the force on the 
top, and is therefore the requested answer. Mathematically, this procedure amounts to replacing M(r) 



X 




M(r) " 
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Fig. 1 
The U-CAS is modeled as a point-like particle of mass m and magnetic moment ^ pointing 

DOWNWARD, HOVERING AT A HEIGHT h ABOVE A NON-UNIFORMLY MAGNETIZED BASE WITH MAGNETIZATION 

DENSITY M(r) PRODUCING A FIELD B(r). 



in Eq.(^ by M(r) = — Mon(6'), where n(6') is a unit vector parallel to n(6'), the latter being defined 
in Eq.(P). The result is 



95, 



dz 



z=h 



X 



— 67rMo desine^/Acos^e + sin^e 
4vr J 

7r/2 
oo 

—dr, 

-h/cos{e) 



for which the r integration is trivial and the 6 integration may be brought to a simpler form by the 
transformation x = sin^(0). This gives |]12 



dB, 



dz 



Ho SttMo 
47r h 



U/A{l-xy + x^dx 



12Mo /jQ 
h 47r' 



which together with Eq.@ shows that 



'11^ 



where 



/o 



477 mg 



is the characteristic length in our problem which will reappear in the next sections. The value of Iq 
may be interpreted as the distance between two colinear dipoles, one of them is of strength /x while 
the other is of strength MqIq, for which the mutual force between them is mg. 

Eq.([TT|) shows that even though the magnetization direction is allowed to vary everywhere inside 



the base, the levitation height h is bounded. Eq.(|ri|) presents the maximum height that can be 
accomplished with a given substance provided that the top is not allowed to move laterally. Clearly, it 
also serve as an upper bound for hm, as stability was not considered yet. Moreover, note that we have 
calculated the maximum magnetic force, which is proportional to dB^/dz. The magnetic field Bz\^^^ 
on the other hand, becomes infinite at z = h, as can be seen by the following simple argument: Since 
each elemental dipole within the plate contributes a field that goes as 1/r^ and since the volume of 
integration goes as r^, the integrand goes as 1/r, for which the integral diverges as r ^ oo. Later we 
show that the divergence of -8^1^=^ implies that the top cannot be stable when placed in such a point. 
In any case, it would be impossible to spin the top. 

C. Uniformly magnetized, infinite base. 

Uniformly magnetized plates are clearly easier to construct. In this section we find out what upper 
bound does this restriction sets on the highest levitation point. In another sense, the result of this 
section may also be considered as a lower bound on the height of levitation (for a top which is guided!), 
provided one considers all possible base configurations. 

We again consider an infinitesimal volume element within the base. It is now oriented along the 
^-direction and interacting with the top's magnetic dipole. If this element is exactly below the top 
(i.e. it is located somewhere on the negative z-axis), it exerts a repelling z-directed force on the top, 
pushing it away. If the element is not exactly below the top then the nature of this force (i.e. weather 
it is repulsive or attractive) is determined by the angle formed between the direction of the two dipoles 
(in this case both point in the ^-direction), and the direction defined by the line joining these dipoles. 
We call this angle 6, and measure it with respect to the positive z-direction, as is shown in Fig.^. 

There exists a critical angle 6q (see Fig.(^), for which the z-directed force vanishes such that for 
7r/2 < 6 < 6q the force becomes attractive. It is therefore useless to put any material within the region 
7r/2 < 6 < 6q and (27r — ^q) < ^ < 37r/2 since this would reduce the magnetic force and decrease the 
levitation height. 

To find 9q we eliminate the z-component of Eq.(^, and find the value of 9 for which it vanishes. We 
pick the value of 6 that lies between 7r/2 and tt. This gives 

^0 = TT - arccos ^/s/S ^ 141". 

We now evaluate the upward force exerted on the top, by summing the contributions of all the 
elements in the base located within 9q < 9 < (27r — 9q). The (uniform) magnetization density within 
this region is taken to be M(r) = Mqz. Using Eq.(0), we find that 



dB,, 



oo 



dz 



= ^lsm9d9 I 2nr^dr\{Moz) 

do -h/cos{e) 

[6 sin 9 cos 99 + 3(3 cos^ 9 - l)r]. 



The integration over r is again trivial, and we are left with an integration over 9. The latter is brought 
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Fig. 2 

For a uniformly magnetized base, all the dipoles that lie within the angular sector 

9o < 9 < {2iT — ^o) (black dipoles), where 9o ~ 141°, repel the top. The dipoles that are outside this 

SECTOR (grey dipoles) attract THE TOP. MATERIAL SHOULD BE PUT only IN THE REPULSIVE SECTOR. 



to a simpler form by changing the 6 variable into x = cos 6, hence 



dB. 



dz 



z=h 



fiQ GtiMq 
Att h 



(5x^ — 3x'^)dx 



h 



-1 

5/2 



This result, together with the equilibrium condition Eq. 
levitation height possible is 



suggests that in this case the maximal 



K 



4.l| 



5/2 



mg Atc 



47r - 



5/2 



/n = 3.5/n- 



(12) 



Note, however, that in order to realize the levitation heights found in this section and in the previous 
one we would need an unlimited supply of magnetic material. In the following sections we find how 
good can we do when the volume of the base is constrained. 

D. Arbitrarily magnetized, finite base. 

Under a given volume Vq, we find the optimum shape and magnetization of the base, that will 
maximize the levitation height. We use the Lagrange's multipliers method to treat this variational 
problem [|l^. 



First, we parametrize the shape of the plate: Let re{6) be the upper integration radius, as is shown in 
Fig-dl)- The lower integration radius is —h/ cos6. Utilizing the cylindrical symmetry of the problem, 
the volume may be written as 



n re{e) 

V = 271 f sinOde I r^dr 

6*1 -h/cose n^\ 

where the angle 6i is determined by the condition that 

reie,) = -h/cosei, (14) 

and is the value of the angle to the upper-right corner of thejkase, as is shown in Fig.(|). 

^T-^ B(r) 




Fig. 3 

A finite BASE IS PARAMETRIZED BY re(6')-THE RADIAL DISTANCE TO THE BOUNDARY FOR A GIVEN 9. INTEGRATION 

IN r IS CARRIED OUT FROM r = —h/ cos9 TO r^ {9)- The polar angle of the upper- right corner of the 

BASE IS 6i. 



The upward force Eq.(|^), which in this case is given by setting 

M(r) = Mon(^), 



10 



is obtained from 



1 47r OB, 



Mq fiQ dz 



■K 

-67r / de sineV4cos4e + sin^e 



(15) 



ex 



cos 9 



re{e) h 

According to the Lagrange's multipliers method |T^, the target function that is to be maximized is 



1 in dB, 



Mq Ho dz 



z=h 



Y, (Vo - V) 



(16) 



where 1/A^ is the Lagrange's multiplier and Vq is the given volume of the base. The target function is 
a functional of rf.{9). We require that a variation of T with respect to Vf^iO) would vanish, thus 



5T 



1 47r (5 



5r,{e) Mofio6r,ie) 
Using Eqs.(|14D,(|15D we find that 

5T 



95, 



dz 



z=h- 



+ 



1 6V 
A^5r,(^) 



0. 



(17) 



sreie) 



o ,,n-nj 3V4cos4^ + sin4^ , r2(^) 
2-71 I d9sm9 < ttttt. ^ ey j 
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X' 



Therefore, the required parametrization re{6) is obtained by equating the integrand in Eq.(|T8|) to zero, 
giving 



X 



3V4:Cos^9 + sm'^e 



1/4 



(19) 



Eq.(^) defines the shape of the optimal base. We see that it is universal in the sense that as Vq is 
changed, the optimal new plate's shape is only a scaled version of the original one. Recall however, 
that 9i is different. 

The value of 9i, determined by Eq.(|T4|), may be combined with Eg . ([T9|) to read 



h 
X 



cos 6*1 



3V4cos4^i + sin^^i = /"^ (^i) 



1/4 



(20) 



which implicitly expresses 6i in terms of h/X. We denote this function by f{x), and write 

Oi = f{h/X). 
In addition, we use Eq. ([l9|) to arrive to an explicit expression for the volume Vq, given in Eq. ([T3|) : 

Vo 



A3 



G{h/X), 



where 



G{x) = — 
^ ' 3 



dO sin 6 



3 V4 cos4 e + sin^ 9 

3 



■13/4 



=/{^0 



V 



+ 



X 



COS' 



^9 
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Similarly, applying Eq.([T9|) to Eq. (|T5|) yields 

dB 



dz 



A = -MoS{h/\) 



(21) 



with the definition 



TT 

S{x) = 6n I desme\/4:Cos'^e + smH 



9i=f{x) 



X 



cos 6 



3\/4 cos^ e + sin^ 6 



1/4 



X 



Using Eq.(^) with the equilibrium condition Eq.(^, we find that 

A 



/n 



Sih/X), 



(22) 



and with few more steps we arrive to 



Vn 



h_hfh 

To~X\X 
h 



'■0 



- = G r; sM ^ 



A 



h 



A 



(23) 



We see that h/lo depends on Vq/I^ through an intermediate variable h/X. It can be solved numerically 
by "running" over a wide range of h/X and evaluating h/lo and Vo/Zg for each of its values. The result 
is given by the dash-dotted line in Fig.(|^). 

Note that at the limit Vq — > oo we find that h = 12/o, in agreement with the result of section( p-B|) . 
The asymptotic behavior oihasVo^O is quite interesting also: This limit may be evaluated from the 
above equations but it is much simpler (and more instructive) to use the following argument: At the 
low volume limit, the base may be considered as a dipole centered at the origin. Such an assumption 
is valid only if {Ii/IqY ^ K)/^o- I^ ^^i^ ^^^^ i^ i^ ^^sy to see that the magnetic force, acting on the 
top, is just 



dB, 



and hence 



dz 



dB. 



z=h 



-/^ 



dz 
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'Att h^ ' 



Att^ h* 



mg. 



Using the definition of Iq we may rewrite the last result as 



We thus conclude that 



h 
To 



lim I — 



.Vo 
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. '■0 J 
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(24) 
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Note that since 



lim ^°/^o 



1/4 



0, 



we see that our assumption is confirmed. The asymptotic hne of Eq. (|2^) is also plotted in Fig. (^ by 
the dash-dot-dotted line. It is conspicuous that it is indeed an asymptotic. 

E. Uniformly magnetized, finite base. 

The solution for this case is essentially similar to the solution presented in the previous section. The 
only difference is in the form of the magnetization that is used. Here we take M(r) = Mqz instead of 
M(r) = Mon(^). The relation between h/lo and Vo/Zg is again given by Eqs.(|23|) with the following 
new definitions: 



G{x) = 



27r 



cos 



1/4 



d9 sin 9 



9 [3cos^(5sin2^-2)] 

' [3cos0(5sin2^-2)]^^^ 



/(^) 



cos^^ 



S{x) 



67r / d9 sin 9 cos 9 (5 sin^ 9-2) 



/(^) 



X 



-1 



[3cos^(5sin2 0-2)] 



1/4 



COS 6' 



X 



The result of this calculation is given by the dotted line in Fig.(|^). Note that in this case h -^ 3.5/o 
as Vq -^ oo, which is again in agreement with the result of section( p-C|) . The asymptotic behavior of 



/i as Vq ^ is identical to the result that was found in the previous section. 

F. Uniformly magnetized, finite, cylindrical base and stability. 

In this case we take the base to be a uniformly magnetized cylinder with a magnetization M(r) = 
Mqz. The radius of the cylinder is R and its thickness is d, with its upper base at the 2; = plane, as 
is shown in Fig.@. 

The magnetic field outside the cylinder is essentially the field outside a solenoid with similar dimen- 
sions. Thus, the ^-component of the magnetic field along the z-axis is given by[|l^: 



^'W - s^^"«' (s4 



(25) 



where 



W {x,y) = 27r 



X 



X + y 



Vx^ + 1 ^y{x + yy + l 
We also define the function F{x,y), to be used later, as 

dW {x,y) 



F{x,y) = 



dx 



We now formulate the conditions for the stability of the spinning top against both axial and lateral 
translations: Stability in the ^-direction is determined by the sign of the 'spring- const ant' in that 
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Fig. 4 
The top hovers stably over a cylindrical base of radius R, thickness d, and vertical magnetization 

Mo- 



direction - k^- Under the adiabatic approximation, the latter is given by the curvature of the effective 
energy along that direction. Hence, 



d'^E, 



eflf 



dz"^ 



z=h 
p=0 



/i 



d^\B\ 



dz^ 



z=h 
p=0 



/^ 



d'B, 



dz^ 



(26) 



z=h 



where p is the radial distance, in cylindrical coordinate system, from the z-axis. 
Similarly, stability in the lateral direction p, is governed by the sign of 



kp — 



d^E, 



eflF 



dp^ 



z=h 
p=0 



/i 



/^ 



AB, 



&" 


B 


' 9p2 


dB, 




dz 


z=h 



z=h 
p=0 



(27) 



1 d^B,. 



2 dz^ 



where in the last equality use has been made of the cylindrical symmetry of the magnetic field and 
the fact that B^ and all of its Cartesian derivatives are Harmonic functions [^ . 

For the spinning top to be stable against translations, both k^ and kp should be positive. Comparing 
Eg. (p^) to Eg. (pG]) , we see that when B^ -^ oo, then kp ^ — kz/2, and therefore one of the pair {kz,kp) 
must be negative. Thus, when the magnetic field diverges at a point, a top placed at that point cannot 
be stable. This proves that the highest hovering height, found in Section ( |11-B| ), is not under stable 
conditions. 

The restriction that both k^ and kp should be positive defines a stable region along the z axis. As 
an example it can be shown |^, that for a base in the shape of a thin disk of radius Rd, the value of k^ 
is positive whenever z > Rd/2, whereas kp is positive for z < ^j2/bRd. The region of stability in this 
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case is therefore -Rd/2 < z < ^/2/5Rd■ Within that region there exists a point Zi for which kz = kp. 
In the case of the disk it is Zi = ^y2/7Rd■ We call this point the isotropically stable point because 
the restoring (stabilizing) force which acts on a top, which is tuned to hover ai h = Zi is isotropic, 
depending only on the deviation from the equilibrium position and not on its direction. 

For a general cylinder we now consider two distinct situations: The first is the one in which the 
stable point is isotropic, i.e. a hovering height h for which kp = k^- The second is the case where h is 
at the verge of stability in the lateral direction. This is also the highest stable point, characterized by 
kp = 0. Using Eqs.(p6D,(p7D we write each of the two distinct conditions as 



-BJ—i, V dz 



dB, 



d'B, 



a 



dz"^ 



(28) 



where a = 6 for the first situation and a = 2 for the second. 

Substituting Eq.(p5|) into Eq.(p8|) defines a functional relationship between h/R and d/R denoted 
by the function G(), 

\^G (d/R) . 

Differentiating Eq . (p5D with respect to z, setting z = h, and using the equilibrium condition Eq.(^, 
gives 

^ = F {h/R; d/R) = N {d/R) , 

where 

N{y)^F{G{y),y). 
Combining it with an expression for the volume, gives 

Va nR^d f d\ /rV d 



vr — — = vr 



-^N^ {d/R) 



Iq /n \Rj \loJ R ' (29) 

G{d/R)N{d/R), 



in R 'i 







which expresses the volume Vq and the levitation height h in terms of a common variable d/R. Running 
over d/R, and evaluating the volume and height according to Eqs.(^), furnishes the required relation 
between h/lo and Vq/Iq. This plot is shown in Fig.(||) where the solid line corresponds to a = 2 (the 
isotropic case) and the dotted line corresponds to a = 6 (the highest stable point). 

Note that both of these plots are not monotonically increasing. They possess a maximum of h at 
some optimal volume Vopt- For a = 6 we find that /imax — l-3/o and Vopt — 310/q, whereas for a = 2 
these are /imax — O.88/0 and Vopt — 90/q. This indicates that using too a much material worsen the 
largest height that can be achieved, which is reminiscent of our conclusion of Sec.( [ll-0|) . If the given 
volume is larger than Vopt however, we can always use only an amount of volume equal to Vopt and 
discard the rest of the material. Hence, in principal at least, one can realize the largest possible height, 
which is why the plots of h vs. Vq had been artificially corrected by assigning the maximum value of 
h for the values of Vq that are larger than Vopt- 
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G. Infinite, arbitrarily magnetized base and stability. 

In this last section we consider the case of an infinite base, which may be arbitrarily magnetized, 
such that the height of levitation is maximized, yet the top is stable against both axial and lateral 
translations. In order to solve this problem, one needs to maximize the magnetic force dB^/dz under 
the constraints that k^ and kp, defined in Eqs.(^6]) and (|27|) respectively, are both positive. The last 
requirement however, results in a non-linear inequality, which cannot be solved analytically. The 
method we take here is to use the constraint kz oc d'^Bz/dz^ = instead, which marks the lower end 
of the stability region along the z-axis. 

Consider an infinitesimal magnetic dipole at the point (r, 6) within the base with magnetization Mq. 
It is situated below the expected hovering position, as is shown in Fig.(^, and its direction makes an 
angle a{r, 9) with the line joining the dipole to the equilibrium position of the top. 




Fig. 5 
The magnetization density at point (r, 6) is Mq and its direction makes an angle a (r, 9) with the unit 

VECTOR f . The angle DISTRIBUTION-a (r, 9) IS chosen so AS TO MAXIMIZE THE LEVITATION HEIGHT WHILE 

MAINTAINING THE TRANSLATIONAL STABILITY OF THE TOP. 



Substituting 



into Eqs.(0),(|), gives 



M = — sin a (r, 9)9 — cos a (r, 9) f , 



dB,, 



dz 



Mo— / d9sin9 

An 



dr- 



2-Kr^ 



(30) 



X 



7r/2 -h/cos{e) 

— 6 sin a sin COS ^ 
— 3(3cos^0 — 1) cos a 
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and 



d^B., 



dz^ 



X 



47r j J r^ 

7r/2 -h/cos{e) 

9 sin 6 (5 cos^ 6 — 1) sin a 

2 



(31) 



— 12cos^(5cos^ — 3) cos a 
The target function to be extremized in this case, is given by 



T[a{r,e)] 



dB., 



dz 



+ A 



d^B. 



z=h 



dz^ 



(32) 



z=h 



and is a functional of a{r,6), with A being the Lagrange's multipher. Since the variation of T with 
respect to a(r, 9) must vanish, we find, on substitution of Eqs.(^) and (0) into Eq. (^2]), that a(r, 9) 
is given by 



2 sin ^ cos + 3- sin 9 (5 cos^ 0-1) 
tan Q;(r,6') = ^ . 

(3 cos2 - 1) + 4- cos 9{5 cos2 9-3) 
r 

Using Eq. (|33|) inside Eq.(^), and requiring that d'^Bz/dz^ = 0, gives the equation for X/h 

TT -008(6*) 



(33) 



d9 sm(9) / ydy 

■k/2 

in which a{r, 9) depends on y according to 



3 sin 9 (5 cos 9 — 1) sin a 
+4 cos 9{5 cos^ 9 — 3) cos a 



0, 



(34) 



.A 



tanafr, 



2 sin 9 cos 9 + 3— y sin 9 (5 cos^ 9 — 1) 
h 

(3 cos2 9-1)+ A-y cos 9(5 cos2 9-3) 
h 



Note that in Eq.(|3^) the variable r has been changed to y = h/r. This way the double integration is 
finite and the singularity in the integrand is eliminated. 
The numerical solution of Eq.(Rl) gives 



A = 0.424/1. 



Using this value inside Eq. (|30D yields 



dB. 



dz 



-1.33 



/io 2'kMq 



4:71 h 



which, together with the equilibrium condition, gives 



K 



8.4- 



4tt mg 



8.4/o. 



(35) 



(36) 



Note that, though the gradient of the field is finite, the field itself diverges, and hence kp = —kz/2. 
As kz = in the case we studied here, we conclude that kp also vanishes. Thus, the lowest stable 
hovering height and the highest one coincide, leaving no range of stability. A small relaxation in the 
conditions however, such as limiting the volume of the material to a finite, though very large value, 
results in a formation of a finite albeit small, range of stability. 
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Base 


Top 


k m 


material 


(5i7)^ax|li 


material 


(BHU^,^ 


p Kg/m^ 




Fe-Nd-B 


320[KJ/m^ 


Fe-Nd-B 


320[KJ/m^ 


7500 


1.35 


Ferrite 


32 KJ/m^ 


Fe-Nd-B 


320 K J/m^ 


7500 


0.43 


Strnatll| 


800[KJ/m^ 


Strnat[|19l 


800KJ/m^ 


8000 


3.2 



TABLE I 

The value of Iq for different base-top material combinations. 



Constraints 


Lower 
bound 


Upper 
bound 


None 
U 

s 

UandS 
V 


1.3/o 
1.3/o 
1.3/o 
1.3/o 
Fig.® 


12/o 
3.5/o 
8.4/o 
3.5/o 
Fig.® 



TABLE II 

Summary of the bounds for the maximum levitation height of the U-CAS under various constraints: 

U=Uniformly magnetized base plate , S^stable hovering, V=under a given Volume 



in. Discussion. 

We showed that the maximum levitation height of the U-CAS is bounded and is given in terms of 
a characteristic length Iq, which depends on the properties of the substance that the base and the top 
are made of. Note also, that n/m = M^^ / p where p is the mass density of the top, and M^^ is its 
magnetization density. Since Mq ~ Br/poj where Br is the residual induction, and since Br is related 
to the energy product [|16| via {BH)^g_^ ~ B^/ApQ, we may also write Iq as 



/n 



po Ml^^Ml"'^ 



{BHrZ.{BH)tl^ 



An 



P9 



TT pg 



(37) 



In this form, /q can be easily estimated from the knowledge of the energy product of the material. 
The best candidates for large hovering heights are the Nd-Fe-B magnets. An example of which is 
Vacuumschmelze's VACODYM 344 HR [0], whose remanence (Br) is 13.5 KGauss and whose density 
is 7.5 gr/cm^. For this magnet, the magnetization is M = Br/An = 1070 emu/cm^. Thus, with 
g = 980 cm/sec^, we find that /o — 1.6 meter. This gives, according to Eq.(^), a maximum stable 
hovering height of about 13 meters!. Typical values of the energy product, density and corresponding 
Iq, for modern commercial magnets available today are listed in Table I. 

According to Eq.(P7D the characteristic length is proportional to the energy product. Therefore, 
the larger the energy product is, the larger Iq will be. A question is then asked: what is the highest 
conceivable energy product? This question has been already discussed by Strnat |jl8|, according to 
which " ...it seems reasonable to assume that the best room-temperature energy products will never 
exceed ~ 100 MGOe", which is about 8 ■ 10^ Joules/m^, and is also included in Table I. These 
predictions, however, assume that a way might be found to give a fairly high H^ to any magnetic 
material. 

Table II summarizes the upper and lower bounds for the maximum hovering height under a selected 
number of constraints. 
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It is important to note, that in the above derivation we assumed that the magnetization vector 
M(r, 6) is independent of the field. Though this is a very good approximation for modern permanent 
magnets hke rare-earth magnets, its only an approximation even for these. Furthermore, the maximum 



field at the top is limited by its mechanical and magnetic strength: Stability considerations show [^ 
that the spinning speed of the top while hovering must be greater than ^y^ Bz\^^j^ /{h ~ h) where 
I3 and Ji are the moments of inertia along the principal and secondary axes, respectively. Also, the 
maximum field at the top is limited by its coercivity as the field is opposite to the magnetization. 
Thus, in practice, Eq.([36|) should be considered as an upper bound for the maximum hovering height 
of the U-CAS. 

Yet another way to increase the hovering height of the top is to use current coil for the base, instead 
of a permanent magnet ^. The advantage of the coil over the magnetic plate is in the fact that one 
can raise the top to the levitation point by electrical means instead of raising the top mechanically, 
as in the permanent magnet case. Here, the hovering height is not limited directly by the strength of 
the magnetic field, but rather by the amount of power that one can deliver into the coil to overcome 
electrical resistance. One might argue that the use of superconducting wires for the coil should lift 
this constraint, but the hovering height is bounded in this case as well, for if the field increases beyond 
a critical value, the superconductor goes into its normal phase. We have shown [^, that for a given 
height of levitation h, the minimum power required for levitation is given by 



A— 



2/^3 



/^m 



Here, P is the power in Watts, r is the resistivity in r2-cm, g is the free-fall acceleration in cm-sec~^, 
/im is the magnetization per unit mass of the top in emu/gr, and A is a number of dimensions 
Amp^emu^/erg^cm^, and is determined by the shape and current distribution of the coil. We found 
that for a rectangular cross-section coil, the minimal value of A ~ 2300. For the optimal coil however, 
we find that A ~ 490. 
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Fig. 6 

Bounds on the maximum levitation height Vs. volume for: Best magnetization possible (dash-dot 

line), Best uniform magnetization (dotted line), Uniformly magnetized cylinder with highest stable 

POINT (dashed line) , UNIFORMLY MAGNETIZED CYLINDER WITH ISOTROPIC STABLE POINT (SOLID LINE) AND 

SMALL VOLUME LIMIT (dASH-DOT-DOT LINE) . ThE TWO CURVES THAT CORRESPOND TO THE UNIFORMLY 

MAGNETIZED CYLINDER HAS BEEN ARTIFICIALLY CORRECTED BY EXTENDING THE HIGHEST LEVITATION POINT TO 

HIGHER VOLUMES. 



